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Br Robert W. Leonard 


SUMMARY 

Williams type modal .solutions of the elementary 
and Timoshenko beam equations are presented for 
the transient response of several uniform beams to 
a general applied load . Example computations are 
shown for a free-free beam subject to various concen- 
trated loads at its center . The discussion includes 
factors influencing the convergence of modal solu- 
tions and factors to be considered in a choice of beam 
theory . Results obtained by two numerical proce- 
dures , the traveling-ware method and JloubolVs 
method, are also presented and discussed. 

INTRODUCTION 

Thu problem of obtaining the response of elastic 
structures to rapidly applied loading is of con- 
tinuing concern to the aircraft industry inasmuch 
as aircraft structures must withstand blasts, land- 
ing impacts, and a variety of other transient loads. 
Tn order to study the various factors involved in 
this problem, it is desirable to consider simplified 
structures for which thorough studies are possible. 
Among the simplest examples of continuous elas- 
tic structures are uniform beams. Consequently, 
beams have been the subject of a considerable 
number of transient response investigations, and 
a variety of solutions of particular beam problems 
are scattered throughout the existing literature. 
(See, for example, refs. 1 to 7. For an extensive 
bibliography, see ref. 7.) 

It is the purpose of the present paper to provide 
a relatively complete source of useful modal solu- 
tions and to discuss the factors influencing the 
covergence of modal solutions and factors involved 
in the choice of the proper beam theory to be used 
in an analysis. To this end, a consistent presen- 
tation is made of Williams type modal solutions 

1 Suporsi'dos \ r A C A Technical No to 4244 by Robert W. Leonard, 1958. 


for the response to a completely general transient 
load of three pertinent uniform beams (a free-free 
beam with a concentrated mass as its center, a 
cantilever beam, and a simply supported beam). 
(Some duplication of the existing literature is in- 
cluded for completeness.) Solutions, based on 
both the elementary and Timoshenko beam theo- 
ries, are obtained by a process which can be 
readily extended to the solution of problems with 
time-dependent boundary conditions. The appli- 
cation of the method is illustrated for the case of 
the free-free beam with a concentrated mass, and 
results for all the beams are summarized in tables 
I and II. In addition, some typical computed 
results are shown for a free-free beam subjected to 
various concentrated loadings. 

Another purpose of the present paper is a critical 
discussion of two numerical procedures, the travel- 
ing-wave method (ref. (>) and Houbolt’s method 
(ref. 8). The procedures are briefly described and 
computations made with both methods are com- 
pared with the modal results. 

SYMBOLS 


effective shear-carrying area of cross section 
arbitrary constant 

propagation velocity of bending discontin- 

IW 

uities, -/ 

\ mr- 

propagation velocity of shear discontinui- 

IA,G 


Young’s modulus of elasticity 

applied concentrated load 

dimensionless applied concentrated load, 

fP 

El 
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( t shear modulus of elasticity 

/ moment of inertia of cross section 
ij integers 

k dimensionless frequency parameter, a?/ 2 

k HI dimensionless rotary inertia parameter, /■// 

k K dimensionless transverse shear parameter, 

i \~kT 

l \ AJi 

I length of beam (half-length in case of free- 

free beam) 

M bending moment (see fig. 1) 

M dimensionless bending moment, Ml /til 

M s dimensionless static bending moment 

nit generalized mass 

m mass per unit length 

m c concentrated mass 

m c ratio of the concentrated mass to total mass 
of the beam, m c jml 
Pi generalized force 

</ applied distributed load (see fig. 1) 

7j dimensionless applied distributed load, 

qP/til 

r cross-sectional radius of gyration 

t time 

transverse shear force (see fig. 1) 

\ T dimensionless transverse shear force, 17 -/ til 

\ dimensionless stat ie t ransverse shear force 
x coordinate along the beam 

y deflection (see fig. 1) 

y dimensionless deflection, y/l 

y { dimensionless translational component of 

/ th natural mode 

y s dimensionless static deflection 

y r dimensionless rigid-body translation 

= V-2 [“ a ' 2 + k «^ + -\l (L 2 -A*«, 2 ) 2 +^- 2 ] 

[«- s 2 +**, 2 ) 

«,*+*,**.* 

7 ‘ p I-Wk* 

5(£) Dirac delta function ^5(f) = 0 for 

r.^-o 

0 dummy variable of integration 



a 

si n & 1 sinh a, 


cos ftd cosh a, 

7z 


£ 

dimensionless space coordinate, xjl 


r 

dimensionless time. 

l 2 \ m 


4>t 

iih generalized coordinate 



rotation of beam cross section 



rotational component of it h natural moue 


static rotation of beam cross section 


O) 

circular frequency of natural vibration 


l(r) 

step function (7(r) — 0 for r<T); 2(r) = 

1 for 


IIV 

© 


Matrix notation: 


[ i 
i i 

I ! 

rectangular matrix 
row matrix 
column matrix 


[ 1 

diagonal matrix 


Prinu 

$ and Roman numeral superscripts are 

used 

to denote partial differentiation with respect 

to 

Dots 

denote partial differentiation with respect 

to r . 




WILLIAMS TYPE MODAL SOLUTIONS 

In normal-mode solutions for the response of 
beams to transient loads, the response is expanded 
in ten ns of a series of normal modes of the beam. 
The coefficients of the expansion (the generalized 
coordinates) are determined from the governing 
differential equations and the boundary and 
initial conditions. Williams type modal solutions 
(ref. 2) differ from ordinary normal-mode solutions 
by virtue of the isolation of that portion of the 
response which may be obtained in closed form by 
a process of direct integration the so-called 
“stati *•” portion of the response. Only the remain- 
ing “( ynamic” portion of the response is expanded 
in series form. 

Tin advantage of the Williams method over 
ordinary modal solutions is its ability to yield, for 
many loading conditions, a more accurate result 
with the same number of terms in the series. 
(See, or example, refs. 4 and 5.) It is particularly 
advantageous where the response function is 
discontinuous. (An example of this is the deter- 
mination of the shear due to a concentrated load.) 
The discontinuity is contained exactly in the 
separated static portion of the response and the 
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series is only required to reproduce a continuous 
remainder. 

In the Williams method, the isolated portion of 
the response is termed static because significant 
parts of the inertia forces are ignored in its 
determination. In general, however, it is time 
dependent by virtue of the time dependence of the 
applied load and of the nonhomogeneous time- 
dependent boundary conditions if such are im- 
posed. In the ease of beams with a fixed point of 
reference, such as cantilever or simply supported 
beams, all inertia forces are ignored in the deter- 
mination of tins static, part of the response; for 
beams with rigid-body freedoms, however, the 
inertia forces due to the rigid-body motion must 
be taken into account. 

One method of obtaining Williams type modal 
solutions is illustrated herein for both t he ele- 
mentary and Timoshenko beam theories. 

EEEM ENTARY BEAM THEORY 

Basic equations. — The motion of a beam 
subjected to an applied load of intensity <[{xj) is 
usually taken to be governed by t he Bernoulli - 
Euler equation 


c> 2 5 2 i/__ 

LI dj 2 ^ m dt i <J 


iu 2 


( 1 ) 


where x is the coordinate along the beam, / is 
time, y(x,t) is the deflection (see fig. 1), EI(x) is 
the bending stiffness of the beam, and m(x) is its 
mass per unit length. The internal bending 
moment M(x,t) and the shear force \ r (x,t) at any 
cross section (see fig. 1) are given by 


and 


M= — El 


d 2 y 

dx 2 




d 

dx 


El 


d 2 y 

dx 2 


( 2 ) 


(3) 


respectively. 

For a uniform beam, these equations may bo 
written in the dimensionless forms 


y lv +y=<i 

(4) 

M=- y" 

(r.) 

V= 

(fi) 



Fhjt're 1. Positive distortions mid positive internal forces 
mid momenta associated with a typical beam element. 


where 


.7 (£> r ) =~'J 
77,. s Ml 


Ji- 


lt I 


V(£j) 


V7 2 

El 


</(£j) 


,,E 

El 


and / is the length of the beam or half-length in 
the ease of a free-free beam. The primes denote 

partial different iat ion with respect to and 

dots denote partial differentiation with respect to 

- t (El 

7 / 2 \ m 

Symmetrical free-free beam with concentrated 
mass. For symmetrical motion of a uniform 
beam having free ends at £= 1 and £= — 1, atten- 
tion is restricted to the portion ()^£^1 with 
boundary conditions stated in the form 


y'( o,r)=o 

(7a) 

y"'((),r)---o 

(71,) 
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»"(1,t)=0 

(7c) 

»'"(!, r) = 0 

(7<1) 


If, in addition, the free-free beam lias a concen- 
trated mass 2m e located at the center £=(), the 
influence of this mass may be introduced into the 
problem by chunking the boundary condition, 
equation (7b), to y'''(0,T) + m f I/(0,T) = 0 where 
m e —m e /ml. On the other hand, the boundary 
conditions, equations (7), may be left unchanged 
and the differential equation (eq. (4)) altered to 

P ir +U +mXZ)]y==q_ (#) 

where 5(£) is the Dirac delta function. In the 
solution that follows, the latter alternative is 
chosen. 

The beam is assumed to be initially at rest and 
undeflected; that is, 

Ptt,0) = Ptt,0) = 0 (9) 

Then the response to a general symmetrical load 
q (£,r) may be obtained in the Williams form by the 
following procedure. 

The solution is assumed in the form 

y(^-r) = ^(r)-|-i/ s (f,r) + Xl^i(r)yi({) (10) 
i=0 

The quantity y r is the rigid-body translation of 
the free-free beam. It is determined to satisfy 
the differential equation 

yr(» = 1 -, 1 =- f q(Z,T)dt (11) 

i -f- m c jo 

and tin* initial conditions 

Vri 0) = »,(0)=0 (12) 

The quantity 2 / s (£,t) is the static deflection de- 
termined to satisfy 

y/' ($,t) = <z(£, 0 _ [l + m c 5(S)]j/ r (T) (13) 

and the cantilever boundary conditions 
y.(0,r)=0 
»/(0,r)=0 
y/'(\,r) = 0 
y/"( i,r)=o 



Note that, by virtue of the definition of y r , y s also 
satisfies the boundary condition y,/"(0,r) =0. 
Finally, tb<‘ shapes i/ f (£) (where i = Q, 1,2, . . .) 
are the natural vibration modes of the beam 
satisfy mg 

Vt tv (i) = V+mMi)\k,*y,(i) 05) 

and the boundary conditions 

Vi( o) = o 
y/"( o)=o 
y/'( D=o 
y’”{ D=o 

where the dimensionless freq uciiey coefficients 

/ jYl 

FI 

Further, it can be shown that the modes ?/,(£) 
satisfy the orthogonality condition 

I 1 [1 +m f «(£)]y # (f)» ; (i)fZ{=0 (tVj) (17) 
J 0 

Note that, by virtue of the arbitrary selection of 
a dat im plane for y s , the dynamic portion of the 
response, in general, still contains a rigid-body 
component (/ — ()). As defined, the total deflection 
y(£,r) satisfies the boundary conditions (eqs. (7)). 
Then remains the problem of determining the co- 
efficients <t>i(r) so that the differential equation 
(eq. '8)) and initial conditions (eqs. (9)) are 
satisfied. 

If expression (10) is substituted into differen- 
tial equation (8) and equations (Kl) and (15) are 
taken into account, the differential equation is 
reduced to 

Td [ii(r)+ki 2 <t>i{T)]y i (^) = — y^ ) T) (18) 

i=0 

Mult plying equation (18) by [1 + m c 5(£)]y^(£) and 
integ ating with respect to £ from 0 to 1 yields, in 
view of equation (17), the following result: 

*,(r) F*,**,(t) = -^ (i= 0, 1,2,...) (19) 

where 



m, = jo + 
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+ w e 6({)]p < ({)y,({,r)rf{ 


beam are self-equilibrating yield, finally, the gen- 
eralized load 


Similarly, substitution of expression (10) into the 
initial conditions (eq. (9)) and taking into account 
equations (12) lead to the following conditions: 


<M0) = 

/\( 0) 
rrii 

(*'=0, 1,2,.. 

-) 

*i( 0) = 

J\m 

m i 

(i=0, 1,2,.. 

► 

.) 


A simple formula for the generalized -mass inte- 
gral Wj for i= 1, 2, . . . lias been presented in 
reference 9 for uniform beams having any of the 
usual end conditions (free, pinned, or clamped) 
but without concentrated masses. In terms of 
the dimensionless quantities defined herein, the 
extension of this formula to beams with a con- 
centrated mass m c at £=() is 

TO f =J^ [l + mVa^)] yH$)d$ 

=4 m c y i 2 (())+ 4 l k 2 l)-2jr/(l)j//"(l) 

+ y t " 2 ( 1)] (i=l ( 2,...) (21) 

For the present case where the end (£— 1) is free, 
equations (21) reduce to 

m i=\ [m c y 1 2 (0)+y, 2 (l)j (?'=1, 2, . . .) (22a) 

The rigid-body generalized mass (z = 0) is seen to be 

w 0 =(l + m f )y 0 2 (22b) 

Some reduction of the generalized -load integral 
Pi{r) may also be accomplislied in general terms 
for ?*— 1, 2, ... . The quantity ll +m c 5(f)]y < (f) 

may be replaced by ^ Vi iv (£) (eq. (15)); then, suc- 
cessive integrations by parts and application of 
the boundary conditions (eqs. (14) and (16)) 
reduce the integral to 

PM = -{* j; V. lv &T)V t (Mi (*=1, 2, ■ • •) 

Substituting from equation (13) and recalling that, 
in natural vibration, the inertia loads on a free-free 


V 2 J„ 2 (*= 1 ’ 2 ’ • • •) ( 23 a) 

On the other hand, for i=0, the quantity Po(t) is 
most simply expressed as 

/ , o(’‘) = 2/oj o y s (t,T)<lz (2:ib) 

It might be noted here that, in the usual method 
of normal modes, the expressions for generalized 
force corresponding to equations (23a) do not 
have the factor 1/Aq 2 . This is one manifestation 
of the more rapid convergence of the Williams 
method. 

The problem now requires direct integration of 
equations (11) and (13) for the deflections y r and 
y Sy solution of equations (19) for the generalized 
coordinates <t> h and solution of equation (15) for 
the natural modes of vibration y u with each 
function satisfying the designated boundary or 
initial conditions. Direct integration of equation 
(1.1) with the initial conditions (eqs. (12)) taken 
into account yields 


Tar. <24) 

Substituting equation (11) into equation (13) and 
integrating four times, taking into account the 
boundary conditions on y s (eqs. (14)), yields the 
following result: 


2/«(£,r) = 


£ ri /'£ 


«tt,r)(rf{) 4 

1 + m e (k~ 6 + 4)J 0 2(t ’ T)rf{ (25) 


The solution of equation (19), satisfying also 
equations (20), is readily obtained by means of the 
Laplace transform. The result is 


) +~ f A(0) sin Ur-e)de 

m i mi j 0 

(i= 0,1,2,...) (26) 

Finally, the natural-mode shapes y t and the 
corresponding frequency equation are derived in 
reference 10. These results, including the natural- 
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f) 

mode shapes and the frequency equation, are 
summarized for easy reference in table 1(a). 

Relations are also given instable- 1(a) for the 
moment A/(£,r) and shear V'(£,r) obtained by 
substitution of the deflection response into equa- 
tions (5) and ((>). (It is also possible to obtain 
these quantities by integrating the total load as 


V=j l (q-y)<ls 



Howi »ver, some care must be exercised in using 
these formulas when the load function is discontin- 
uous in time or has discontinuous first derivatives 
with respect to time.) 

Other configurations. -The response of a 
uniform free-free beam without a concentrated 
mass is given by the results in table 1(a) with 
7 r c ~(). The response of a cantilever beam may 
also be obtained from the response of the free-free 
beam with the concentrated mass by a limiting 
process in which the mass ln c approaches infinity. 
Results for the cantilever beam are summarized in 
table 1(b). For completeness, the Williams 
solution for a simply supported beam is shown in 
table 1(c). 

Time-dependent boundary conditions. It is 
worthwhile to point out that the method outlined 
in this report is directly applicable to the solution 
of problems with nonhomogeneous time-dependent 
boundary conditions. Such problems require the 
separation of the solution into two parts, one 
satisfying the time-dependent boundary conditions 
and the other capable of being expanded in terms 
of time-independent functions such as the natural 
inodes of the beam. (See, for example, ref. 11.) 
In the Williams method, this separation is already 
made and time-dependent boundary displacements 
or forces are simply introduced into the boundary 
conditions imposed on T/ s or into the equations for 
rigid-body displacements 

Consider, for example, a uniform beam fixed at 
one end and given a variable displacement at the 
other, such that its differential equation and 
boundary conditions are 

y IV (£,T) + ylt t T) = q{t,T) 
T/(0,r)=--p / (0 l r) = y"(l l r) = 0 


?/( 1 ,t) = < 7 (t) 

The solution would be assumed in the form of 
equation (10) but with y r = 0 since there is no 
rigid-bod}^ translational freedom in this case. The 
static portion of the solution would be determined 
to satisfy 

yJ v (£,T) = 7 i(£,T) 

and 

y,( 0 ,r) = y/ ( 0 ,r) = y" ( 1 , t) = 0 
y,( 1 ,t)=//(t) 

while the expansion functions y { (where i= 1,2,.. .) 
are the solutions of 

vt IV (t)=kt*yi(i) 

and 

2/t(0) = Pi / (0) = y<"(l) = i?i(l) = 0 

(the natural modes of a clamped-pinned beam). 
In order to complete the solution, the generalized 
coordinates corresponding to a beam initially at 
rest and unstressed would have the usual form 

<Mt)= - />,(r -+ f r /’-(8) sin a-,(t -e)de 

m, m , J u 

(i— 1,2, . . .) 

where 

/'.(t-) 2 | q{S,T)y,(Z)il$ 

w,=J o y ?(£)</$=— 2 ]/ 2 y/0)y/"0) 

Similarly, a uniform free-free beam with a speci- 
fied tirre-dcperident displacement </(r) at its renter 
moves iccording to 

2? zir (£,T) + 2/(f,r)=5(f,T) 
y( 0 ,r) =y (t) 

y'( 0 f r) = y n (\ y T) = y f "(\ t r) = 0 
In this case 

Thir) 

and ?/,(£, t) is determined from 

T ) '/(£’ r )“(/(0 

and 

^0),r) = y/(O l T) = y "(l > T) = y/''(l,r)-0 
wliile the natural modes ly, (where /=!, 2 , . . .) 
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TABLE I.— RESPONSE OF A UNIFORM ELEMENTARY BEAM TO A GENERAL LOAD 


(a) Symmetrical free-free beam with a concentrated mass 


Quantity 

V(S,t) 

Analytical expression 

y r (T)+yAS,T)+^Z 4 >iW)yA& 
1 = 0 

M(£,t) 

A/ s ( £ , T ) + <t>i(T) Mi ( £ ) 

i= 1 

V(£,t) 

1 r Ak* t) + S 4 >i( t) Vt (£) 

i = 1 

4 >i(r) 

i 

- /,i(r )+— ( T /W) sin £,(t-0) <10 

m, m t Jo 

I'M 

' 

fc 1 * Jo* (i=U 2 ,.. .) 

f/o f yAS,r)di (>=(>) 

do 

\ 1 
j 

! Mi 

^ [wi f yj(0)+y, 2 (l)] 0=1, 2, . . .) 

(H-mJy„ 2 (i= 0) 

T)M 

. ' — ( [ f </(£,r)f/{(</T) ! 

1 + ^h Jo Jd Jo 

yAS,r) 


mam 

1 -J'j;' 5 («,r)(rf*)>+ r+ L(f-{+ 0 J>(i.'»''f 

j V.(i,r) \ - f*j({,T)rf{+ fi- 1 ) g(«,r)rf« 

J 1 l “T /w r J() 

i 

HAS) 

[cosh > cos %/<’,£+ cos \ fri cosh \ k& 

Ific 1 / • , • , ; w > , 1 , ^ /sin \fr,£— sinh \ /f,-£ , cosh a A*,£— cos \ k^\] 

— (sm >A-,+smh (cos \ A* J --|-c(>s!i A J 1 - . ... . . , + , , , T 

- \ sin \ A^+sinh \ k t cos ^■A’f-reosn yA*, / J 

MAS) 

— ( r ki [—cosh \k f cos \ &t£+cos A*, cosh \ A*,£ 

, m 0 , • , , w , , , 7 \ /sin N A\-£+siuh \ A*,£ cos \ jfc^+eosh a A’,£\~| 

+ - 0 \A^ (sin^A^ + sinhyArJ (cos N A%A<-osh> A%J 1 . . — ---------- 

^ \ sm \ A’i + sinh ^ k ( cos -v -cosli \ k f / J 

e,(f) 

— f r Av V2 [cosh sin v A\£+cos \ A, sinh \ A*,£ 

| , 777 r ,7-/. , , ,-w , , , iy\ /cos v t ( £ + (*osli A A^£ sin > A’ ( £— sinh \ 

! + t \ki (sm \ A",-p smh \ kj) (cos^ /r, + cosh \ k,) 1 . _ - . - T - + , - , , , I 

^ \ sin > A*,-|-sinh k,- cos ^ A^ + cosn \ k, / J 


Frequency equal ion: 

cos A k sinh A fc+sin \k cosh \ k-\-m v yk ( 1 -f- cos \ k cosh \ k) ---0 


; 50 — 2 


50202 :' 
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TABLE I. — RESPONSE OK 


A UNIFORM ELEMENTARY BEAM TO A GENERAL LOAD — Continued 


(b) Cantilever beam 


Quantity 

V&T) 

) 

PM 

Thitr) 

M s (£, r) 
r,«,r) 


Analytical expression 


y «(& t) +S^f( r )2/(© 

;' = l 


+ ^>2<t>,(T)M t (£) 

i = l 


V.(tr)+^(r)V t {& 


1 f,' ( C((0) sin k t (r— 6)<ld 
in i ui i j 0 


a ! 2 




U) 


n r* 


0 Jo J 1 J 1 




- J { «(f,r) (r/£) 2 

— J q(t,r)<iz 


2b(£) 

^ / sill -y A-,£— sinli vAr,*- 
\ sin -y/ti+sinh \A", 


/ sin V*il+ sinli W ,£ 


\ sin yAy + sinh y7 * 

r,(£) 

^ ■/. :i/2[ rosVArif+cosh -\ A-,£ 


\ sin y Ay + sinli y iy 


cosh VAy£— cos \ Ay£ 
cos yfly + eosh yAy 


cos y Ay£-f cosh yi.f 
cos yAy + cosh y/y 


sin y Ay£ — sinli \Ay£ 
cos VAy-f- cosh \ Ay 


1 + cos \ A' cosh yA:— 0 


Frequency equation: 
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TABLE I.— RESPONSE OF A UNIFORM ELEMENTARY BEAM TO A GENERAL LOAD— Concluded 

(c) Simply supported beam 


Quantity 

Analytical expression 

yd,*) 

<M r ) sin 

i= I 


+S i 2 7T 2 0, (t) sill 7 7T £ 
1 = 1 

mr) 

V t s ($,t)+S iVViW eos 

i=l i 

<t>, (r) 

— 2/h(r) + 2 i 2 7 r 2 | I\{0) sin i 2 7 r 2 (r— 6)dQ 

Jo 

P i ( 7 ) 

If 1 . 

•4 4 Sill |7t£ f/f 

'hr 4 J 0 

Ih (£,t) 

f £ f * f £ f £ 5(£, 7 ) w?) 4 -i (TIT w© 4 +5 (1 ~* 2) IT ^ ,r) (,/y2 

Jo Jo Jo J o Jo Jo Jo Jo 0 J o J o 


I - f f 9 («,r) (</*)’+$ \ ' f * ®(€,r) CrfO* 

Jo Jo Jo Jo 

V,%t) 

— f 9(f, 7 )<^+ | 1 9 (i.O (rf€)* 

Jo Jo Jo 

i 


are the modes of a cantilever beam. Or if, instead 
of y(0,r) —g(r) } there is given the force boundary 
condition y //f (0,r) = A(r), the rigid-body motion 
is determined from 

Vt{t )= ( q{^r)dl— A(r) 

Jo 

The static solution y s is taken to satisfy 

y/ v {t,T) = q(Z,T) — y r (r) 

^*(0,r) = y/(0.T) =p/ / (l ,t) = 2//"(1,t) = 0 
and the modes l/< (where r=0, 1,2, . . .) are the 


natural modes of a free-free beam. In this case, 
it can be shown by integrating the differential 
equation governing y s that y/ ,/ (0,r)=A(r). 

Thus, the treatment of problems with time- 
dependent boundary conditions involves no spe- 
cial separate procedure when the Williams* method 
is used. 

TIMOSHENKO’S BEAM THEORY 

Basic equations. — In the elementary beam 
theory, deflection occurs only by virtue of the 
rotation of the beam elements and only their 
translational inertia is taken into account. The 
Timoshenko beam theory (ref. 9) permits addi- 
tional deflection due to transverse shear and 
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accounts also for the rotational inertia of the beam 
elements. According to this theory, the motion 
of a beam subjected to an applied load of intensity 
q(x,t) is governed by the equations (see, for ex- 
ample, ref. 6): 



where ^ is the rotation of the cross section (see 
fig. 1), r is the radius of gyration of the cross sec- 
tion, and A S G is the stiffness in transverse shear. 
The effective shear-carrying area A s differs from 
the total area because the shear stress is not con- 
stant over tlie cross section. The bending mo- 
ment M and transverse shear force V are given by 


Symmetrical free-free beam with concentrated 
mass. — For the application of Timoshenko's 
theory to the symmetrical motion of a uniform 
free-free beam with a mass 2 m c at the center £=0, 
attention will again be restricted to the portion 
1. As in the case of the elementary theory, 
the effect of the mass may be introduced into the 
differential equations if desired. However, for 
illustrative purposes, the differential equations 
(30) will be left unchanged and the mass will be 
introduced in the boundary conditions; the 
boundary conditions then become 

*( 0 , t )=0 

F0,r)— iK1,t) = 0 
yi j/'(0,r) = m r 2/(0 , t) 



M— — EI (28 ) 

Of 

and 

v - a -° (&-*) < 2,j > 

For uniform beams, these equations may be 
written in the dimensionless forms 

2 (y' ~ 

I. liv' -4')' 

A, s 

(:n) 





when' the transverse shear coefficient 

is a measure of the freedom of the beam to deflect 
through transverse shearing action and the rotary 

inertia eoeflieient kiu — j is a measure of the rota- 
tional inertia per unit length. 


Note that the functions y and are both neces- 
sary for adequate definition of the deformation of 
the beam. Since these go hand-in-hand, the terms 
“solution" and "response," as used herein, will 
apply to these functions collectively and the single 
notation 2/(f,r) ; ^(£,r) will he used to specify both 
functions. 


Note that t he location of the concent rated mass 
at £ = 0 and the restriction to symmetrical motion 
excludi any effect of the rotational inertia of the 
eon cen rated mass. 

The beam is assumed initially at rest and unde- 
flected: hence, 

y(t o)=(T 

y(i, o)=o 

► (34) 

*(*,0) = 0 

km= o. 

With the problem thus completely defined by the 
differential equations (eqs. (30)), the boundary 
conditions (eqs. (33)), and the initial conditions 
(eqs. (14)), the solution may be obtained as 
follows 

Assume that 

Ff r ) = F(0+F($>0+S$<(0yt(£) 

i=o 

where p y (r) is again the rigid-body translation of 
the beam, y«(£,r) ;^q(£,r) is the static solution, 
and T/i |);tfq(£) (where i — 0, 1,2,...) are the 
natural vibration modes. Tin* rigid-body trans- 
lation of the beam y r is governed again by the 
differential equation 
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and the initial conditions 


y r (0) = y r (0) = 0 (37) 


The static solution is determined to satisfy 


t/'+i k <$;-*.)= o 

Kg 

jri (y*-4>,Y=-~i+y’ 

ft s 


A 


(38) 


equation (43b) by k RJ 2 \f/ J} adding the two equations 
and integrating the sum over the range 
yield the result 


&+*,**= (i=0, 1, 2, . . .) (44) 

Hi j 

which takes into account the orthogonality rela- 
tion (eq. (42)). The generalized mass and general- 
ized load appearing in equation (44) are 


and the cantilever boundary conditions 

y>{ o,r)=o 

*,(0,r) = 0 

1^/(1, t) = 0 

y/( l,r) — tA*( 1 >r) = 0 - 


(39) 


The mode shapes y ( (£); ^ ( (£) satisfy 
**"+%* 

7T2 (Vt’ —ftY +k f *y t =0 

the boundary conditions 

*i(0)=0 

*/(i)=0 

2//(l)— ^ ( (1)=0 

tt 2 lit' (0)=— m e k,*y i(Q) 

H-s j 

and the orthogonality relation 


(40) 


(41) 


( * {[l + iMtf)]^)^) 

“b&*/Vi(£)^/(£) }d£=0 (i 9^j) (42) 

The derivation of this orthogonality relation is 
shown in the appendix along with the solutions 
to equations (40) and (41). 

Substituting equations (35) into the differential 
equations (eqs. (30)) and utilizing equations (38) 
and (40) reduce the differential equations to 


co 


XI (4>i+ kr<t>,)y i= — y s 

i = 0 

(42a) 

X (0i+*(V«)^i=— s 

i =0 

(42b) 


Multiplying equation (43a) by [1 +m f 5 (£)]*/., and 


and 


*i (j)=j ^ [l + m c 5(f)]?/iy s +A! ft rVi^^^ 


respectively. 

By a similar process the initial conditions (eqs. 
(34)) become 


MO) 


MO) 


Pj( 0) 

Pi(0) 

Mi 


(i=0, 1,2, . . .) I 


(i=0, 1,2, . . .) J 


(45) 


It is shown in the appendix that the generalized 
masses m t of the given free-free beam can be 
evaluated for z=l, 2, . . .as follows: 



(i=l,2, . . .) (46a) 


The remaining generalized mass m 0 reduces, as in 
the elementary theory, to 

w 0 =(l+ i c )f/(i 2 (46b) 


since the symmetric rigid-body mode has no 
rotational component ^ 0 . Further, the general- 
ized force 1\(t) for i 1 , 2, . . . may be reduced, 
by a process of substitution from equations (38) 
and (40) and integration by parts, to 


P i (r)= /Jo 


(i= 1,2, . . .) (47a) 


and the rigid -body generalized force P 0 is 


r 0 (r) = 3/" yJt,T)(it 


(47b) 


On the basis of the assumed form of the solution 
expressed in equations (35), the problem of deter- 
mining 9/(£,t) ; ^(£,r) has been replaced by a number 
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of component problems requiring determination 
of the functions //,+), yj.£,r) <k(r), *uid 

ViiV ;«/+.£)• The solutions of these component 
problems must now he obtained. 

For j/ r (r), integration of equation (30) in con- 
junction with equations (37) yields 


'/,(>) = , U- | | J «({,r)f/{ («/r)* (48) 

1 -f //ifjo J<) Jo 

For y, £,t) ;^(f,r), substituting ?/ r from equation 
(30) into equations (38) and integrating, in con- 
junction with the boundary conditions (eqs. (39)), 
gives 


V*{t,T)= ( * f * [ * Vq(£,T) (dZY-k* ( 5 | € 5(f,r)(r/f)* 

Jo Jo Jl Jl Jo Jl 

+1 +»,[)'• (I’-O-dr !+?)] £ 




-xtj: 


«G.r) (^) 3 - 


rk<s-H)x> 


r)</{ 


49) 


Since equations (44) and (45) are identical to 
equations (19) and (20) of the elementary solu- 
tion, the generalized coordinates </>,(t) are again 
given by equation (20). The solution of equations 
(40) for the natural modes ydX) + t (£) is given in 
the appendix. 

The solution just obtained and corresponding 
solutions for cantilever and simply supported 
beams are summarized in table II. 

TWO NUMERICAL METHODS OF SOLUTION 
TRAVELING- WAVE METHOD 

A traveling-wave method for calculating the 
response of a structure to transient loads cun be 
devised only if the motion of the structure is 
governed by differential equations of the hyper- 
bolic* type. The simplest beam theory which 
completely fulfills this requirement, is Timo- 
shenko’s theory, which includes the effects of both 
transverse shear and rotary inertia. 

In developing a traveling-wave method of solu- 
tion, it is convenient to first replace the Timo- 
shenko partial differential equations (eqs. (27)) 
with the following four equivalent ordinary 
differential equations written along four charac- 


teristic lines I + , I — , 
/, t plane: 

Along 1+ where 

As Ci 


II + , and II 


— AM-\- mr 2 dSl— Vdt 

Ci 


the 


0 

(50a) 


Along I— where * 


—AM — m r 2 A U + V A t — 0 

Cl Ci 

(50 b) 


Along 11+ where ( ~=~- 
as c 2 


AV 


Along II- 


— mdr-r (mc 2 U-\-q)At — 0 (50c) 


1 


when* ( f-= — — : - (IV 
(Is C 2 C >2 


+ m As d - (m c 2 12 — q) At = 0 (50d) 

The derivation of equations (50) is given in refer- 
ence 0. The dependent variables are the moment 
M, the shear V, and the linear and angular veloci- 


ties v~- 


and c. 


0/ aMd 


The quantities 


FI 


Iaj ; 

- V" >n 11 


\ m r 2 

are the propagation velocities of 


discon inuilies in moment and shear, respectively 
(phase velocities of disturbances with infinitesimal 
wavelmgth). In each equation, the total differ- 
entials specify infinitesimal differences in the 
design. ited characteristic direction. 

For any given beam, the slopes of the charac- 
terislit lines are known; hence, closely spaced 
netwoiks of eharaeteristic lines may he drawn in 
the sp. ice-lime plane. Various schemes for the 
approximate step-by-step integration of equations 
(50) o' or such networks arc possible. In general, 
all require some form of interpolation since 
Timoshenko’s equations have two characteristic 
nets. (The particular case where the two nets 
coined e. c A — c 2 , has been treated in detail for 
uniform beams in ref. (j.) One integration scheme 
is briefly described in this section. Attention is 
restricted to a uniform beam for which the char- 
acteristics are straight lines. 



ON SOLUTIONS FOR THE TRANSIENT RESPONSE OF BEAMS 


13 


TABLE II.— RESPONSE OF A UNIFORM TIMOSHENKO BEAM TO A GENERAL LOAD 
(a) Symmetrical free-free beam with a concentrated mass 


Quantity 

A 1 1 aly t i cal e xp ressio 1 1 

y(i,r) 

2/r (*r) +y*(f ) 'r) +S <Mt) */*(£) 

j — 0 



m, T ) 

MAS,t)+±, <j>i(T)M t (£) 
! = 1 

V(Z,t) 

C.(£,t)+S </>, (t) r,(£) 
1 = 1 

4>, (t) 

- / - ,? (r)- + — ( T /’iW Sin k,(T—0)de 
mi m t Jo 

/Vt) 

| ' h^rJhAm (i— 1 , 2 , . . . ) 

f/„ ( yAt,T)<lt (i= 0) 

Jo 

j 

VI i 

(l+m r )y 0 2 C- 0) 

~Jr (t) 


V.&r) 




■V x ((,r) 


j 

t ****** «-■>/>« 

Vi(& 

V. A “I 

-I, ^( a . 2 +pi)~' ( cosl ‘ a ‘+^ cos A) [j'os A£— cosh a^+A, (sin P£—J ( ^ sinh 
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TABLE IT.— RESPONSE OF A UNIFORM TIMOSHENKO BEAM TO A GENERAL LOAD— Continued 
(a) Symmetrical free-free beam with a concentrated mass — Concluded 


Quantity 

Analytical expression 

JTr 

s/Pi 2 — W\ f Pi sin Pi . , sinlia, . 

( { K 7f a Sinh oc,$+ s in 

V 0i / ^ a f 0 t oti 

(8 2 —k 2 k 2 ) r / b \~ n 

2 J ( ros * 1 <*i+yi ( ‘ os &) 1 A<(cos Pit— cosh a<$) — ^sin 7i sinhc^J 1 > 


-(’W-kft*) { 7, 0i cosh a ^ + SU ^ COS P A 

h Pi Oti 

+ m c ^ ( ( ' os, > “<+"r< cosPi ) Qv.^sin Pi^+^ sinh a f ^+ (cos P,^+y t cosh 

M,(H) 

r t (& 

r ,kr f PiSmPi . . sinli . 

0 4 sinh a,£-|- sin 

Pi «( Pi «< 

— m, ^ * (cosha,+Y< cos/3,) j~A,^cos /3 J+^ cosli «,£^— ^sin Pi(—-*‘ sinh a t $yj\ 

A, 

a 

i sin Pi sinh a t 

oti 

cos Pi~\~ - cosh a, 
7/ 

7 1 

Pi'-W* 

i 

«i 

9 (£% 2 + A^/ 2 ) -f- -y (A:* 2 A*/?/ 2 ) 2 H“^, 2 J 

Pi 

A-y 2 L^+^+V *» 2 - fc ^ 2 ) 2 +/t 4 J 

Frequency equation: 

7 sin (3 cosh a+cos p sinli a+m f ^ ^ ) sin 0 sinh a+^7+^ <*os 0 cosh aj— 0 



ON SOLUTIONS FOR THE TRANSIENT RESPONSE OF BEAMS 


15 


TABLE II. -RESPONSE OF A UNIFORM TIMOSHENKO BEAM TO A GENERAL LOA !)■ —Continual 

(h) (-ant Hover b(»ani 


Quantity 

Analytical expression 

j/ft.O 

/=l 

*(f,0 

f= 1 

Pft.r) 

TfJi,T)+it<t>i(T)Mi(£) 
i = \ 

Off, r) \ 

(= i 

4>i (V) 

— / ^ T ' + — ( /L{0) sin kfir — d)<l6 

ni f mj „ 

Piir) 

^i| o ' 

m j 

| (1 ' |pl*ftH-*«Y<*ft)] 'ft=~ - 2 J. *f(l) [JJ,. *'(1)1 4 

VAt,r) 

{ { { { ) £ j; «ft,r)(rft) 4 -*,* j £ J’ £ 5({,r) (r/£) 2 

'Ps(^t) 

( ( ( 8 ft, T) (</{)* 

-O s (f,r) 

- rj £ 7^.0 ('ft )2 

T.ft.r) 

-J 8ft,T)fft 

y,(& 

('[As j8,{— cosh a,i+ A, ftn ~ A sinli 

+*(& 

a* a* 2 A*, 2 r / \"i 

^ ^ |^A/(cos ft£ — cosli of/{)— ^sin 7 f sinli «*£ jj 

Mi(S) 

^ (ft- 2 ki 2 h% 2 ) j^Af ^sin fijt ;- sinh a^^-b(POs ft£+7/ eoslt 

T,ft) ! 

( 0- [ A - ( < OS < ' os *' “I?) - (sin P£— ^ s >nli "'(=)J 

A, 

Pi ■ i 

Sill ft — - Sill ll ai 

a, : 

<■ os ft H cosh atj 

y ; 

7* 

aP+kW 

ft*-* f **,*_ 

1 

<Xf \ 

ktyj I [- <*.*+*«/*) + yC 2 -/.-^ 2 ) 2 +p] 

Pi 

2' [(*,*+**/*) + ^W-*«, 2 ) 2 +p] 

frequency equation: 

2+^ sin P sinli a “b^7+ ^ cos p cosh a— 0 


50292 ; 


59 — 
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I A 10, K II. — HKSPONSK OK A I'XIKOR M TIMOSIIKXKO UK AM TO A OKXERA1, LOAD 

(<■) Simply supported l.eain 


Quantity 


T7 (£,r) 

V(S'T) 

1'At) 


IfJZ.r) 


Anal ylieal expression 



//,«($. C) +S |0/ (t) If .,,, - 

r l0i(r)lt. = ».y- sin /V£ 

<A(£d 

l "f(l — ^' 2^2 ) [0/(r)]j- r 

+ (d ^f)[0 ( (r)],, 


.\/»({.T) + X) )'V -^(l— 1— ^^2 )l0i(r)U.-/,.^ > 

!’*(£, t)+XJ d -^«i S f0>WJ*..,, l( +6 i 2 [0((T)]t., A j- cos nr£ 


sin nr£ 


1 ‘M.k, r„ 


m j m 


/> 


0) sin kt(r — 0)de 


1 r> _ 

2 I </0;,t)<u\ /tt£ (f£ 


,£ r r r r r r r r p 7 <v> ^ 

U Jll Jo Jo J(> Jo Jo Jo Jo Jo 




*.<' f.r) I 1 | S I *5(«.r)(r/t)»- | ' | * ( ? | * tfft.r) (</£) '-(jf- 1 ; ) | ' ( £ 7 (f.ri <//{) 2 

Jo «/o Jo Jo ./i) Jo Jo \- n/ Jo Jo 


MaS.t) 

T’AZ.t) 

Frequency t'cjua lions: 


£ v($.t) (,/£)*- i $ ( ( * 7 U.O O/?) 2 
0 */ 0 J 0 J 0 


| £ <l(Z,T)<i$+ | 1 ( £ 7(i,r)(^) 2 
Jo J) Jo 


A’i = </, or />, where 

//, — ^ ^ 1 d v 2 7T (/%“ + A'«/") — \ [ 1 ~h £V 2 ( i"* 2 -)-/'/*/ 2 ) ] 2 — 4i* N 2 A*«^ 2, f 4 7r 4 ^ 

- ^ -y/;*- >^1 + iWOcA+kmA + \ |1 4- tV 2 (t s 2 +A- fi/ 2 )] 2 — 4A* s 2 A-« ; ? / V^- 
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TIip boam is divided arbitrarily into segments 
Xr (see fig. 2), and the time interval is taken 

according to A/=— Xr. This specifics a lattice of 
c-i 

points in the 1 space-time plane at tin* intersection 
of the' 1+ and I— characteristic, lines. Consider 
a general point 1 (fig. 2) from which character- 
istics of both families have been drawn backwards 
in time. The 11+ and II— lines have steeper 
slopes (since o 2 <+i) and terminate at points 2' 
and 4'. Then the differentials in equations (50) 
may be replaced by the appropriate finite differ- 
ences and tlu‘ following equations are obtained: 

1 ( M, - 1 /.,) 4 m r- (12, — S2 a ) — A * (1',-M ',) - 0 (5 1 n ) 

<’i - ' 2 

1 ( A/, - .1 /.,) - mr 2 (12, — S2 4 ) + % ( C, 4 l \) = 0 (5 1 1 >) 

c i - 

1 (Ti— IV) — mO’i — r-y) 

C>> 

-f- me, A y (a, + 12,.) 4+ (f /1 -T </.') =- 0 (51 c) 

- (1 i l 4'J "4" m 0’ i /’i 1 ) 

C 2 

- me*-* (12,-1 Uv)- A ' (</.4 f/v)- OJ (51(1) 


It is assumed that 1/, \\ i\ and 12 are known 
at points 2, d, and 4. Parabolic interpolation 
formulas are substituted into equations (51 c) and 
(51 d) to give the quantities at points 2' and 4' 
in terms of their values at points 2, 3, and 4. 
Then ('({nations (51) become four equations for 
the four unknown quantities ;\f u 1+ r,, and 12 t 
in terms of known values of M y \\ r, and 12 at 
each of the lattice points 2, + and 4. These may 
be solved to obtain a matrix recurrence formula; 
however, a simplification may be introduced based 
on the fact that quantities at points 2 and 4 have 
already been determined to satisfy the. following 



x 

FifiCRK 2. (irkl scheme for traveling-wave numerical 
procedure. 


(‘({nations: 

- - (Af«-Af<) - mrHil ,- <>,) + A ' (1',+ (’.) =() 

c i ' ' - 

1 (Af t - a/*) + »»•* (i 2 4 — - + 1 + r.,) = o 


(52) 


(Combining (‘({nations (52) with (‘({nations (51a) 
and (51b) hauls to the result 

mr-mu- — 0',- r,) --A) (55) 


where M has b(‘en eliminated. Now (‘({nations 
(51 e), (51d), and (5A) const it ill e thri*(‘ (‘({nations 
for the three unknowns !2 t , \\, and r, in terms of 
known values of 12, \ \ and r at each of the points 2, 
3, 4, and 5. Solution of these equations yields the 


recurrence formula 



r 

u. 


a* 

;«4 



7i 

A 

1 

1 • A' 2 " 

kU 

V, 

+i*tj 

v 3 

- i*ei r« 

+i.i,i 

T\+\q] 

72 

73 




r o 





t' r, 

74 

> 


( 54 ) 
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<*+$)%* 


_d<p 

dr 


and 

2 K 

C-, 

K 

‘i 

Cl 

k ri 

C-? 

r, 2 

1 

C 2 

Cl 2 

kjii 


v 'lien 


„ 2 


^-kuKd+K*) -lt KI (\+K*) %(\ 4 -IO) 






-K* 


K 

<'\ 

0 


1 

0 


0 

0 

11 K\ 


[-UI • 


(*-*$) 


V er /Cl 


£i K 

Cl 


±K % 


c 2 


Ci K 
<‘i km 

cl 1_ 

Ci* kjti 


1 k RfK (1 A 2 ) k M ( 1 -j" A 2 ) ■ - 2 ( 1 + K 2 ) 

Cl c 1 



-1 

— - K 

Co 

0 

[.'Id- 

C >K 
c 1 

K 2 

0 


0 

0 

0_ 


k R! K 2 


I 


k'lfiK 1 


IQ\- 


; y k R! -^K 
y c i 


J k(\ +k 2 


; y k M * c -*K(\+K*) 

~ r, 


1 / /' 
«- /i>/ A 

- f’i 


( 1 -J*I) WA'iH A’ 1 ) J *«* ~ 2 A’ (1+A 2 ) 


l_c, A/ 

k„*c t 2 


Note that the more appropriate parameter ^ ^ ias usot l here instead of k„. 

The response of a beam may now be obtained by the repeated application of equation (54) except that, 
as is indicated in figure 2, special formulas which take into account the particular boundary and initial 
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conditions of the problem, must be derived for 
boundary and initial points. In addition, it must 
be remembered that the characteristic lines are 
possible loci of discontinuities in the dependent 
variables or their derivatives. (See, for example, 
ref. 6.) Such discontinuities will arise, for 
example, if there are concentrated loads (or imposed 
velocities) that have histories which are discon- 
tinuous or have discontinuous derivatives. Dis- 
continuities in M or 12 are propagated with veloc- 
ity Ci (the locus of such a discontinuity is shown 
by the dashed characteristic lines in fig. 2) ; discon- 
tinuities in V and v propagate with velocity c 2 . 
The magnitudes of discontinuities may be prede- 
termined through the application of equations (50) 
in a manner which is illustrated in reference 6. 
Thus, in general, discontinuities may be added as 
they are encountered in the step-by-step solution. 
In the scheme which led to equation (54), it will 
be noted (fig. 2) that special consideration is 
necessary for points just above the locus of a dis- 
continuity. Special formulas are certainly re- 
quired where there is a jump in one of the functions 
and it may be desirable to account also for dis- 
continuous first derivatives. 

Once 12, V, and v have been determined at a 
point, M may be obtained at that point by in- 
tegrating equation (50a) or equation (50b) along 
the proper characteristic from some boundary 
whore M is known. 

It should be mentioned that the selection of the 
1+ and I — lines as the basic network is based on 
(rather intuitive) considerations of the stability 
of the numerical procedure with regard to propa- 
gation of errors. (See ref. 12.) It is assumed that 
the domain of dependence imposed by a numerical 
procedure should at least encompass the total 
theoretical domain of dependence. This would 
not be the case if the steeper 11+ and II — char- 
acteristic lines were utilized as a basic grid. 

HOUBOLT’S METHOD 

A cantilever beam acted upon by a series of 
concentrated loads p u j > 2 , . . . }>„, is shown in the 
following sketch. 

P\ P 2 P ? * ’ * P n 


* L ■ Ax -) 

0 I 2 3 ' ’ ’ n 

Station 


Such a beam has the deflection 

l?/r| = [tf]l2>l (55) 

where the subscript c is used here to indicate that 
the deflections are cantilever deflections measured 
with respect to station 0 and where [f?] is a matrix 
of stiffness influence coefficients. The inverse 
equation is 

bl =[£]-%<! (56) 

The deflections of a free-free beam may be ex- 
pressed in terms of the cantilever influence coeffi- 
cients [(?]. For the free-free beam the symmetric 
deflection is y—yo+Ve where y 0 is the deflection at 
station zero (the center of the free-free beam); 
hence, equation (56) becomes 

1 

1 

\p\ = [G\- 1 \v\-vJLG\- 1 ; (57) 

1 

But, on the free-free beam, there is the additional 
load po . From the condition of overall equilibrium, 
p 0 is given by 

i = 1 

i: 

i 

+]/oli i . * . UK?]- 1 ; (58) 

i 

and equations (57) and (58) may be combined 
into a single matrix equation for the loads at sta- 
tions 0 to n of a free-free beam. This equation 
may be written 

l/>l=MM (59) 

where 





Ill 

[b\ = -[l 1 . . . 1J[6T 
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1 

1 


a = [\ 1 . . . UK?]- 1 • 


whore 


1 


\n = lH)l ml 


and whore the vectors ]/j| and \y\ now contain p {) 
and //o terms, respectively. 

Consider the symmetric motion of a free-free 
beam subjected to the applied distributed load 
(/UJ). The equivalent concentrated loads at the 
stations 0, 1, ... n are, at any time, 

l/>l = ar|(|f/| — Imllyl) (60) 

where [ni] is the diagonal mass matrix, and [\V] is 
the weighting matrix (ref. 1^)* 


Thus, die Houholt method simultaneously deter- 
mines all the deflections y(Xi,tj) (where i = 0, 1, 
. . . //) in terms of the deflections y{x u tj^x), 
y(Xi,tj - : .), and y(x h tj^). A concentrated load 
2f(t) al the center of the free-free beam may be 
included in equation (02) by adding to the rigid- 
hand side the term 

w\ m 

whore 


7 0—1 

2 20 2 


2 20 2 


2 20 2 
2 20 2 
— 1 6 7 

The esseidial feature of the Houholt procedure 
(ref. S) is the method of expressing t he second time 
derivative y. The acceleration at time tj=j&t is 
writ t en as follows: 


1 

0 

0 


The -solution of a problem by repeated applica- 
tion of equation (02) requires that the initial con- 
ditions of the problem at time / = 0 be expressed 
in tern s of fictitious ordinates [?/|_ 1 ,|?/|_ 2 at times 
t. ] and This is accomplished by expressing 
the time derivatives | ?/’ 0 and |i/| 0 in terms of a third- 
degree curve passed through poirds at t— t u 0, 
t-i, and t„ 2 at every station. For the case where 
the beam is initially at rest and tin* applied loads 
increase continuously from zero initial values, the 
initial < onditions 


\<>\ [ mi + (A % 


l?/l«=|y|o=|j/:o=« 


(64) 


'-^•'“(ao* (( >d 

and is obtained by passing a third-degree curve 
through the points at f—tj, tj- U tj ~ 2 , and tj-i at 
each station jq. Substituting equations (bO) and 
(bl) into equation (59) (written for time f — tj) 
and solving for leads to the recurrence equation 

l.vi>=[^]k/lj+ ( ^ )2 [C)(5|i/k- l -4|j/k_ 2 +|i/|,_3) 

(62) 


are converted by this process to 
!?/ ! o — 0 ' 

?/!-,=- — ?/! V 
.?/!-*=— 8 1 /!,^ 


( 65 ) 


Then, equation (62), applied for 7=1, yields 


\y\i=[BY<i\i- 


M* 
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which may ho solved for |?/j, to obtain 


'[ |,|+ (i? |r| ] |B|1 


The application of equation (62) for /— 2, II, . . . 
is now straightforward. All that remains is the 
determination of the cantilever influence coeffi- 
cients [G], 

In this connection it should be pointed out that, 
although no restrictions have been made on the 
beam theory to be used, the application of the 
method, as formulated, with Timoshenko’s theory 
requires that the deflections ?/ be interpreted in a 
general sense and include also the rotations of the 
cross sections \p. Thus, with two quantities to be 
determined at each station, the order of the 
matrices is 2 n and the required computational la- 
bor is roughly four times that required with the 
use of the elementary beam theory. A compro- 
mise which affords increased accuracy over the 
elementary theory, yet avoids this large increase 1 
in computational labor, is the use of a theory 
which contains transverse shear freedom but no 
rotary inertia. With no associated inertia load- 
ing, the rotations \p need not be explicitly included 
in tlie step-by-step dynamic analysis and do not 
appear in the recurrence formula, equation (62). 
The cantilever influence coefficients are deter- 
mined as follows on the basis of this latter theory. 

The influence function (Green’s function) 
G(x;ji) is the solution y of the equations 




[ l/; (2 -*)] — 


and the boundary conditions 


y(o)=^(o)=^(/)=s:(0-^(/)=o 


where G /t and G s are the solutions of the differen- 
tial equations (equivalent to eqs. (67)) 


d 2 , d 2 

5/ 2 ^ di 2 * fl ) 


l A -° 


and boundary conditions (equivalent to eqs. (68)) 


Gn (0;-r,)= ' (0;/,)=^ 

_ 0 . v d 2 6V, ,_ n 

d.r LI i>x 2 1 0 

G s ( d ^(/;,r,)=0 


For given distributions of bending and shear stiff- 
nesses, equations (70) may be integrated directly 
in conjunction with the boundary conditions (eqs. 
(71)). For a uniform beam, the resulting total 
influence functions arc: 


i/; ‘ Ei\ 2 o / >,r 


In dimensionless terms, equations (72) may be 
writ ten 


<?(«;«.) (i.-0 (*<*.) I 

=*.*«!+ 2* («“|) 


r (78) 


and, for a uniform beam, the recurrence formula, 
equation (62), becomes 


However, the deflection may be written as the 
sum of bending and shear contributions, ?/=?/*+?/, s , 

with and it is expedient to write for the 

influence function 


|/ >’ | -[l->l i |'H ] 


G (x;j\ ) = G a (x;jci ) + G s (x;x t ) 
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1 

1 

5=li i . . - lpr ; 

1 

1 

1 

|6|=-[^r i ; 

1 


DJ=-U 1 ■ • • U [OV 1 



The elements G i} of the matrix of influence co- 
efficients [6 T ] are seen to be 




(Kj) 

=k*j&i + 1 ( 

;-0(A£) 3 

(>>j) 

‘signates the 

row and j, tin 

s column. 


RESULTS AND DISCUSSION 

MODAL SOLUTIONS 

For illustrative purposes, example computations 
have been made for the case of a uniform free-free 
beam, for which A*/e/=0.1, £ s =0.2, and m c ={). 
The beam is subjected to the applied concentrated 
load 

5(£,t)==2/(t)5(£) 

wher e/(r) has each of the time variations shown in 


figure 3 The given values of k RI and k , are appro- 
priate to a beam having a solid rectangular cross 
section and a ratio of half-length l to depth equal 
to 2.887. The calculated response has been 
limited to the history of transverse shear at the 

point £=-. Moment calculations arc omitted 
because they do not provide as severe a test of the 
analytical methods. The point is chosen 

arbitrarily, since the location corresponding to 
maximum transverse shear is not known in 
advance. 

The espouses to the step and ramp-platform 
loads were obtained on the basis of elementary 
theory from table 1(a) and on the basis of Timo- 
shenko’-* theory from table 11(a). In each case, 
six modes were used in the expansions. The 
responses to the other three functions (figs. 3(c), 
3(d), aid 3(e)) were obtained by superposition of 
the step and ramp-platform results. The re- 
sulting shear histories are shown in figures 4, 5, 
h, 7, aid 8 up to a time corresponding approxi- 
mately to the period of the first natural mode of 
vibralicn of the beam. 



0 I .2 .3 0 .1 .2 .3 



0 .1 2 .3 0 .! .2 3 



T 

(a) Stop function. (b) Ramp-platform function, 

(c) Square pulse. (cl) Triangular pulse. 

(c) Mast pulse. 

Figure IL — Some fundamental load histories. 
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The long (lushed curve in each figure is the 
static portion of the response. Thus, the largest 
dynamic overshoot factor, achieved by the step 
loading (fig. 4), is approximately 2.7. 

It will be noted that the blast pulse load (fig. 
3(e)) has only one-half as much impulse as that 
contained in the square and triangular pulse 
loads (figs. 3(c ) and 3(d)). Hence, if the responses 
to the three pulse loads (see figs. 6, 7, and 8) are 
to l)e compared on tin 1 basis of equal input impulse, 
the response to the blast pulse must be doubled. 
On this basis, the blast pulse is seen to cause- the 

highest shear stress at 

An indication of the convergence of the modal 
results in figures 4 to 8 is provided bv the bar 
graphs in figure 9. On each graph, the heights of 
the bars correspond to the magnitude of the static 
portion of the response (zero frequency) and to 
the amplitudes of the terms in the series expansion 
for the dynamic portion. (For each load, r is 
sufficiently large so that the load function has 


attained its constant value.) The bars, thus, 
represent the maximum possible contribution of 
each term to the total. Note that, for the loads 
of long duration (the step and ramp-platform 
functions), the static part and the first term con- 
tribute a proportionately large share of the re- 
sponse and sufficient accuracy could he obtained 
with only three modes. On the other hand, the 
adverse effect of reducing the load duration is 
illustrated by the bar graphs for the responses to 
to the three pulse loads. In each case, no con- 
vergence is apparent for the* first few modes and, 
in the eases of the square and blast pulses, there 
is some doubt as to the adequacy of the even six 
modes, particularly with the use of the elementary 
theory. 

Further evidence, of the effect of load duration 
on convergence and, in addition, an indication of 
the effect of load distribution are given by the 
following eases of the response of a uniform simply 
supported beam to various loads: 



24 


TECHNICAL REPORT R-21 — NATIONAL 


-< cn 


AUTICS AND SPACE ADMINISTRATION 


Elementary theory (6 modes) 

Timoshenko theory (6 modes) 

Traveling- wave method (6 stations) 



free beam to a ramp-platform load concentrated at £ — 0. 


Case (3): For a uniformly distributed impulse 
load q(£ t ) = 5 (t), 

\1 (£,r) =A \ sill i.7rf sill iVr 

K *=1,3, . . . t 

\ r (£,r) —4 X] cos iirf sin tVr 

i= 1,3. . . . 

Case (4): For an impulse load concentrated at 
t lit* cenliM- q(£,r) = d ($— ^ 5(r), 

i — 1 

M £,r) — 2 X) (—1) sin ?*7 t£ sin iVr 

1 = 1,3, . . . 

t-1 

5!j ( — 1) ? cos iTrf sill iV 2 r 

i = 1, 3, . . . 

These results are based on elementary beam theory. 
(The last two cases are simply time derivatives of 
the first two cases.) The effect of load duration is 
illustrated by comparing cases (1) and (S) (or 
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cases (2) and (4)). Changing from a load of 
infinite duration (the step load) to a load of zero 
duration (the impulse) introduces a factor i 2 and, 
hence, reduces the rate of convergence (and, in 
fact, produces divergent series in both case (3) 
and case (4)). Similarly, the effect of spatial 
distribution may be seen by comparing eases (1) 
and (2) (or cases (3) and (4)). Changing from a 
distributed to a concentrated load introduces a 
factor i and hence reduces the rate of convergence. 
The apparent change in the sign of half the terms 
is not significant since each series is essentially an 
irregularly alternating series (except at certain 
specific combinations of £ and r). 

THE NEED FOR TIMOSHENKO'S THEORY 

The question of which theory should be list'd to 
determine the response of a beam to a transient 
load is intimately related to the convergence of 
the result. This is because tin* secondary effects 
of transverse shear and rotary inertia become 
increasingly important for the higher modes. 


Note in figure 9, for example, the growing disparity 
between the natural frequencies of a uniform free- 
free beam calculated on the basis of the elementary 
and Timoshenko theories. Thus, if it is deter- 
mined that, for a given beam subjected to a certain 
load, modes strongly affected by transverse shear 
and rotary inertia contribute a large share of the 
response, it is unlikely that elementary theory will 
yield correct results. In case the given beam is 
a complicated non uniform structure, a rational 
procedure for determining the proper theory would 
be to consider a uniform approximation to the 
given beam, quickly obtain the response of the 
uniform beam to the given load by elementary 
theory (table I) and investigate the convergence of 
the response, and at the same time to consider the 
influence of rotary inertia and transverse shear on 
the modes (as manifested by the differences in 
natural frequencies obtained with the elementary 
and Timoshenko theories). 

This reasoning is generally, though not- conclu- 
sively, confirmed by the results in figures 4 to 9. 
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Only the first mode of the free-free beam shows 
good agreement between the frequencies as given 
by the elementary and Timoshenko theories. (See 
fig. 9.) Thus, the responses to the stop and ramp- 
platform loads (figs. 4 and 5), which depend 
heavily on the first mode and the static contribu- 
tion, also show good agreement between the two 
theories. The responses to the pulse loads (figs. 
(>, 7, and 8), obtained with the two different 
theories, bear little or no resemblance to each other 
since they depend heavily on the higher modes. 
However, except for the response to the square 
pulse load (fig. (>), the two theories do yield about 
the same peak stress. (It is felt that the positive 
peak achieved in the first half period should bo 
given more weight than the negative peak achieved 
later, since the latter would be considerably di- 
minished by material damping which has not been 
included in this analysis.) 

Boundary conditions also influence the need for 
a more refined theory. For example, although the 


elementary theory is adequate for obtaining the 
shear < ue to a step load on the free-free beam of 
figure 4, it cannot be used to obtain the shear at 
the cei' ter of the same beam where the input is a 
prescri )ed ‘‘step-velocity” of the 1 point £— 0 (the 
so-called problem of the “instantaneous arrest of 
the root, of a moving cantilever beam”). In the 
latter case, elementary theory yields a divergent 
series (ref. 14), whereas Timoshenko’s theory yields 
a modi 1 solution which converges to finite values 
(ref. <>). 

NUMERICAL SOLUTIONS 

The two numerical procedures, the traveling- 
wave method and ITouboll’s method, have also 
been used to calculate the transient response of 
the un : form free-free beam considered in the pre- 
vious sections to an applied ramp platform load. 
In the calculations by both procedures, the beam 
was divided into six segments ( 0.1 (>(>7). Re- 
sults are shown in figure 5. 
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For the traveling-wave method, the time inter- 
val is 


At = £*,A{ = 0.01<W7 


Relatively simple boundary formulas (based on 
linear interpolation) were used in this ease and 
discontinuities in the derivatives of 1 and r, arising 
from t he discontinuities in the slope of the ramp- 
platform function, were ignored. 

in contrast with the traveling-wave method, the 
Houbolt method imposes no inherent restriction 
on the selection of the time interval At in relation 
to the space interval A£. The time interval may 
be taken as large as is consistent with the desired 
accuracy. This freedom has been utilized in that 
t lie calculations by the Houbolt method have been 
made with the time interval Ar = 0.(W!W3, which is 
twice the time interval used with the traveling- 
wave method. In addition, the computations by 
the Houbolt method do not include the effects of 
rotary inertia and were stopped at a point just 
beyond the peak load. 


It will be noted that tin* numerical results in 
figure f> approximate the (essentially converged) 
Timoshenko modal solution fairly well. Both 
numerical methods underestimate the peak stress, 
the traveling- wave method by 14 percent, and the 
Houbolt method by 9 percent. 

The greater accuracy of the Houbolt result is 
particularly significant since it was obtained with 
less computational labor due to the use of the 
larger time interval. It must be kept in mind, 
however, that economical use of the Houbolt 
method requires that it be applied in connection 
with the elementary theory or with the addition 
of transverse shear alone. If both rotary inertia, 
and transverse shear must be included, tin* labor 
required in applying tin* Houbolt method is quad- 
rupled. Fortunately, rotary inertia is negligible 
in many problems. (See, for example, ref. l. r >.) 

It should hi* mentioned that the traveling-wave 
procedure, as so far conceived, has proved some- 
what sensitive to minor changes in tin* scheme 
used to obtain recurrence formulas. For example, 
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Elementary theory Timoshenko theory 

Ay = 0 5.59 30.2 74.6 139 223 326 0 5.03 19.9 36.3 50.6 58.8 68.6 
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the simplest and most obvious scheme (in which 
the 11+ and II— characteristics are extended 
backward from the point 1 (fig. 2) to intermediate 
points on the I— and 1+ characteristics connect- 
ing points 5 and 2 and points 5 and 4 and linear 
interpolation is used to determine the unknowns 
at these intermediate points) lias been found to 
yield significantly less accurate results. A pos- 
sible reason for this sensitivity stems from the fact 
that the slope of the basic network of characteristic 
lines is dependent on rotary inertia. In fact, 
rotary inertia is necessary to give the beam equa- 
tions the wave character essential in the concep- 
tion of a traveling-wave method. Thus, in view 
of the relative negligibility of rotary inertia for 
many practical problems, this sensitivity is per- 
haps not surprising. In general, it must be con- 
cluded that a traveling-wave numerical method of 
analysis which is superior to the Houbolt method 
has not yet been devised. 

Numerical results have also been obtained for 
the response to t he triangular pulse load by super- 
position of the ramp-platform results. These 
results are shown in figure 7. The results indicate 
that more degrees of freedom must bo taken with 
both procedures to predict adequately the response 
of the beam to the given triangular pulse load, and 
from this example it appears that the modal 
method of solution is to be preferred. However, 
the simple problems discussed herein do not por- 
tray the main advantages of numerical methods. 
For example, numerical methods are readily ex- 
tended to apply to nonuniform beams and con- 
veniently adapted to the use of modern high-speed 
computers. A fundamental characteristic of nu- 
merical methods is the replacement of professional 
engineering time by routine computing time. 


Hence, numerical procedures are not to be con- 
demned on the basis of the results in figure 7. 
The selection of the best method requires the con- 
sideration of all these factors in relation to the 
specific problem. 

CONCLUDING REMARKS 

Williams type modal solutions, based on both 
the elementary and Timoshenko beam theories, 
have been given for the response of several uni- 
form beams to a general transient load. The 
response to any specific load may be obtained 
from these solutions by performing a series of 
indicated direct integrations of the load function. 
Typical computed results have been shown for 
the shear response of a free-free beam to various 
concentrated loads. 

The convergence of modal solutions is shown to 
depend both on the history and distribution of tin 1 
load. Decreasing either the duration of the load- 
ing or the region over which the load is applied 
reduces the rate of convergence and may produce 
divergence. 

The need for a more refined theory, as compared 
to elementary theory, is intimately related to the 
rate of convergence of the modal solution. If 
modes which are strongly dependent on trans- 
verse shear and rotary inertia contribute a large 
portion of the response, Timoshenko's theory must 
be used. 

Comparison of the Houbolt and traveling-wave 
numerical methods indicates that the Houbolt 
procedure has many advantages over the traveling- 
wave procedure as so far conceived. 

Langley Research (’enter, 

National Aeronautics and Space Administration, 
Langley Field, Va., February 5, Wotf. 



APPENDIX 


SYMMETRICAL NATURAL VIBRATION OF A UNIFORM FREE-FREE 
TIMOSHENKO BEAM WITH A CONCENTRATED MASS 


NATURAL MODES AND FREQUENCIES 

Tile differential equations and boundary condi- 
tions governing symmetrical natural vibration of 
a free-free beam with a concentrated mass at its 
center may be written in the following dimension- 
less forms: 


2 ( y * id ± k 2 k — 0 

a 

- (Al) 

-<f)'+k‘ 2 y=() 


J 


H o)=o 

(A2a) 

*'0 )--=<> 

(A 21 >) 

37' ( l)-*(l )--=<> 

(A2c) 

y' (0) + 7 n f fc*y( 0) =0 

(A2d) 


where //(£) ; ^(£) is the natural mode and k is pro- 
portional to the circular frequency of vibration. 
Each of t lie solutions of equations (Al) has the 
form 

y(Z)^A<* 


X= ±?’0 where 

:> [ (*.*+**/*) +\/ (L 2 -***)^ £] 

r 

:i [fe’+y)+-y/ (*,*-*,/*)*+£] 

(A 5 ) 

The general solution of equations (Al) may be 
written in the form 

y (£) — ( ' cosh «£ + f 2 sinh <*£ 4 ( ’ 3 cos 0£+ f \ sin 0£" 

^2 I 7,2 L 2 

*(0 = 1 - (<% cosh af +<A ftinh a£) 

a " V 

-K - * * (rw-os 0f-r 3 sin 00 

(A8) 

where f , f 2, ^3, and an* arbitrary functions of 
A: which must be determined in order to satisfy 
the boundary conditions (eqs. ( A2) ) . 

Subst tuting equations (AO) into equations 

(A2) yields four homogeneous algebraic equations 

for (\ 2y f’ 3 , and (\. The existence of non- 
trivial solutions of these equations requires the 
vanishir g of the determinant of coefficients. This 
criterior yields the frequency (‘({nation 


where A and B are arbitrary constants. Substi- 
tuting this form into etjuations (Al) leads to the 
biquadratic (‘({nation 

X 4 ±i*“(A’A-j k'/ir 2 )\* -1 k 2 (k 2 k s 2 kf i 2 — 1)— 0 (Al$) 


and to the following relationship between A and 
B: 


B 



(A 4 ) 


Equation (AB) has the four solutions X=±a and 



7 sin {3 cosh a + cos j 3 sinh a 


+ 


where 


(a 2 -\ k 2 k. 2 ) f / a 0\ . . , 

a'-yfP L- H (ra.) S, ' ,/S8ml, “ 

+( 7 +^) cos /3 eosli al=0 (A 7) 

_a 2 +k*k* 

7 ~P-Ie 2 lc} 


The solutions of equation (A7) are the natural 
frequencies of vibration if, where f=0, 1 , 2, . . . . 
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For k—k u the homogeneous algebraic equations may be solved for the relative magnitudes of the 
quantities C u C 2t (\, and (\, For i— 1 ? 2, . . the resulting vibration mode shapes may be 
written in the form 


_ , , ... f sin Bt , si nil a t 

= ( A r~ cosh a,Z— ----- COS 
L Pi 01 i 


f-t 2 b 27. 2 r V 

TO,' ' * (cosh a, + 7 , cos ft) cos A{-coah «,( 

<*i +/?i L 

+A, ^sin siiih 0= 1,2, . . .) 

, /lA kfk* (fr sin/?,: . , „ . s 

+t(Q = ( R A „ y ‘ a — smh a,£+ 

P, Pt 


Sltlll a i . — Pf—kfkJ 2 , , 

sin 2 [o*i ( ( * os ‘ 1 a i 

oci aS-f-ps 


+7 * cos /SdjjMeos /3,£— eosli «*£)— (sin 7* sinh 


( sh 


«,«)]} 


0 = 1,2, . . .) 


y (a 8 ) 


where 


8. . . 

sin fit sum a* 

€t i 

cos Pi 4 — cosh a f 
y t 


The rigid body mode, corresponding to fco=0, 
has the components 7/ 0 (f) = t r and ^<>(f)= 0. 

ORTHOGONALITY OF THE NATURAL MODES 

The differential equations (Al) are satisfied by 
any of the infinite number of natural modes and 
corresponding frequencies. Tims, for the ith mode, 

*,”+ £ 2 (y/-*.)+*«tt«Vi= o 


TJiis process is valid also if the roles of the i th 
and jth modes are reversed. Interchanging i and 
j in equation (A9) and subtracting the result from 
equation (A9) hauls to 

(kf— k?) Jo (T/iT/j l ^‘1 ( Th , —'h)Vs 

~h (A10) 


(?// ^i)' J rkt 2 yi — 0 

Let the first of these equations be multiplied by 
the rotational component \ pj of the jth mode and 
the second by the translational component y } . If 
the resulting equations are added and integrated 
over the beam length, there results 

A-, 2 J (j?,y,+ J o (?7i' — h)'VA 

- f h’^A— \ it* (Vt'-hWA 

Integrating by parts the first two integrals on the 
right-liaml side of the preceding equation yields: 


If now the boundary conditions (cqs. (A2)) an* 
imposed, equation (AlO) is found to reduce to 


Jo 


(AV 


or 


A'/)TO ( t/((0)y i (0) 



[1 + m r 5 (i)]j/i(?) 3 /^(i) 
+A’*/Vi(e^j(e J- ^{=0 


since 


OVA) (All) 


J mM£)yt(Z)yM)<lS= 


to f y, (0)j/,(0) 


fri 2 J (i f/ i y i +k B i l 4'fl'j)di= — [j. 1 2 (y/—h) ?/>+iA , .V,J n 
+ f (A 9 ) 


Equation (All) is the orthogonality condition 
satisfied by the natural vibration modes of a uni- 
form free-free Timoshenko beam with a concen- 
trated mass at the center. 
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DETERMINATION OF THE GENERALIZED MASS 

The determination of the generalized mass 

m i == j +^r5(£)]5i 2 (£)+^ie/Vi 2 (S)^rfS 

by direct integration is a somewhat laborious 
process for 0. Fortunately, m* can be ex- 

pressed in terms of certain boundary values of the 
mode shapes by the application of a limiting 
process to equation (A10) in which the functions 
\li and \pi are considered as continuous functions 
of k. Thus, if 

fc j— lc i~\- dk 

and 

yj=Tti+>iy=yi+^<ik 

4 , j=4 , i i +^j- 

equation (A 10) becomes, in the limit as dk ap- 
proaches zero, 

J>, [£, (5/-W (!), 



(A 12) 


This equation is applicable to uniform beams and 
may be extended, if desired, to nonuniform beams. 


On substitution of the boundary conditions (eqs 
(A2)), equation (Al2) reduces to 



+Mm=-m c ym 


(i=l,2, . . .) 


Hence, he generalized mass is given by 

(A13) 

Note that only the second boundary condition 
(eq. (A2b)) is altered by differentiation with re- 
spect to k. This arises from the fact that only 
this boundary condition depends for its satisfac- 
tion on the frequency equation (A7); that is, the 
mode shapes (eqs. (A8)) satisfy the other boundary 
conditions for any value of k but satisfy equation 
(A2b) only for k = k t since 

*'(i)=-V-* 2 LW) 

a 


where F{k) is defined by equation (A7). It is 
only by virtue of this dependence of one or more 
boundary conditions on the frequency equation 
that equations (Al2) yield a value of /«,. Hence, 
it must be concluded that, for a beam for which 
none of the boundary conditions depend for their 
satisfaction on the frequency equation (as, for 
example a simplv supported beam), equations 
(A 12) are not applicable. For such a beam, m t 
is determined by direct integration. 
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